Measurement-induced two-qubit entanglement in a bad cavity: 
Fundamental and practical considerations 
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An entanglement-generating protocol is described for two qubits coupled to a cavity field in the 
bad-cavity limit. By measuring the amplitude of a field transmitted through the cavity, an entangled 
spin-singlet state can be established probabilistically. Both fundamental limitations and practical 
measurement schemes are discussed, and the influence of dissipative processes and inhomogeneities 
in the qubits are analyzed. The measurement-based protocol provides criteria for selecting states 
with an infidelity scaling linearly with the qubit-decoherence rate. 
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I. INTRODUCTION 

Entanglement is one of the key features of quantum 
mechanics, and during the past decades it has been 
demonstrated experimentally in many different physical 
systems. By coherent control of interacting quantum sys- 
tems, entangled states can be engineered directly [lHaj - 
Alternatively, entanglement can be established as a con- 
sequence of the outcome of a measurement process — ei- 
ther as a continuous (possibly quantum-non-demolition 
(QND)) measurement [3 or as a consequence of a sin- 
gle quantum jump Some of the above examples 
employ cavity quantum electrodynamics (QED) [10( for 
mediating the interaction between the quantum systems, 
which allows for the direct engineering of entangled states 
in the strong-coupling regime [l|, [H, @| ■ The present pa- 
per considers a different case — the bad-cavity limit of 
QED, in which the damping rate of the cavity field is 
fast compared to the coupling rate between qubits and 
the cavity held. Hence, any information of qubit co- 
herence being encoded into the cavity held will immedi- 
ately be lost from the cavity and the above-mentioned 
direct-engineering schemes are inapplicable. However, 
turning to a measurement-based protocol, the detection 
of a held transmitted through the cavity allows for re- 
establishing a hrm knowledge of the qubit state and hence 
for the creation of entangled states through measure- 
ment back action. The measurement is of the continuous 
type, which is theoretically well-described by stochastic 
master-equation methods [llj ■ In contrast to the work of 
Ref. [l2j using similar theoretical methods, our calcula- 
tions are not restricted to the dispersive and linear regime 
of the coupling between the cavity held and the qubits 
but allow instead for a more generalized set of parame- 
ters (even a resonant coupling) in search for the optimal 
choice for entanglement generation. We consider feasible 
experimental approaches and discuss the physical limita- 
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tions imposed both fundamentally by the measurement 
process and practically by decoherence mechanisms. 

This paper is arranged as follows: The basic idea 
for the protocol is outlined in Sec. [TTJ while the the- 
oretical modeling is elaborated on in Sec. IIII1 Various 
practical measurement schemes are analyzed in Sec. IIV1 
while qubit decoherence and inhomogeneities are added 
in Sec. [Vj After a general discussion in Sec. I VII we sum- 
marize the conclusions of the paper in Sec. I VIII Some 
mathematical details are deferred to the appendix. 



II. ENTANGLEMENT-GENERATING 
PROTOCOLS: THE BASIC IDEA 

The physical setup under consideration (see Fig. [1]) 
consists of two qubits placed in a cavity subjected to an 
external driving held and to a continuous measurement 
by employing a phase-sensitive detection of the held leak- 
ing from the cavity. Let the ground and excited states of 
either qubit be denoted by |g) and |e), respectively, and 
consider the two-qubit basis set (Jee) , |gg) , |+) , |— )}, 
where |±) = (|eg) ± |ge)). Our aim is to generate 

the spin-singlet state, |— ), by a probabilistic detection 
scheme with a high fidelity and a high success probabil- 
ity. This state does not couple to the cavity held when 
the coupling parameter g is equal for the two qubits. 
The other states, |ee), |gg), and |+) span the spin-triplet 
space and the cavity held may induce rotations within 
this subset of Hilbert space. Furthermore, the coupling 
between the qubits and the cavity held gives rise to a 
correlated decay mechanism, which induces transitions 
(with rate 27 P ) within the triplet space: |ee) — > |+) and 
|+) — > |gg), while the singlet state, |— ), is unaffected. 

The idea is now to prepare a separable initial state, 
and subsequently to collapse (probabilistically) \ip) into 
|— ) by the measurement process. The separable initial 
states of opposite spins, |eg) = ^(|+) + |— )) or |(|g) + 

|e»(|g> - |e» = I |gg> - I |ee) + ± |-), both have a 
50% overlap with the desired singlet state, and our task 
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FIG. 1. The physical system under consideration. Two qubits 
are coupled to a cavity field, a c , which is driven externally by 
a constant coherent field, ft (in the frame rotating at the driv- 
ing frequency, ujl)- The field-decay rates through the cavity 
mirrors are denoted by m and K2. The output field, At, is 
subjected to a balanced homodyne measurement using a lo- 
cal oscillator, qlo = |oiLo|e lS . The differential and summed 
photo-currents are denoted by i-(t) and i+(t), respectively. 



is to identify a measurement scheme, which is able to 
distinguish between the singlet and triplet components 
of and thus to facilitate the state collapse. 

No real or virtual transitions can take place within the 
one dimensional singlet subspace of the qubits and it ex- 
periences no interaction with the cavity field. The triplet 
subspace, however, consists of three states, and their in- 
teraction with the cavity field can be tailored to affect 
the transmitted radiation field in two different ways: the 
cavity field and the damping of the spin system via the 
cavity mode can drive the collective spin to a steady state 
mean spin polarization which causes a phase shift of the 
transmitted radiation, or the dynamical driving of the 
triplet spin components can induce a frequency modu- 
lation of the transmitted field auto-correlation function. 
The steady state change in the transmitted field is visi- 
ble in the homodyne photo-current, while the frequency 
modulation can be observed with a lock-in detector. 



III. THE STOCHASTIC-MASTER-EQUATION 
APPROACH 



Our knowledge of the physical system depicted in 
Fig. [T] is accounted for by the density matrix, p, which 
evolves according to the stochastic master equation 
EE!: 



J2 V 



[c m ]p + ^jn[d}p((t), (i) 



where the Hamiltonian, H , describes the interaction be- 
tween the cavity field, the qubits, and the coherent driv- 
ing field, /3j n = /3e~ l " L *. Decay processes are modeled 
by the super-operator, V[c m ]p = - \& m Cmp- \p& m c m + 
Cmpc^, while our knowledge from the continuous moni- 
toring of the system is incorporated by the measurement 
super-operator, T-L[d]p = dp+pd) — (d+dft) p. The detector 
quantum efficiency is denoted by 77, and the real- valued 



function, £(t), models the randomness of the detection 
process with ensemble characteristics, (£(i))E = and 
(£(*)£(^'))e = S(t — t'). In the frame rotating at the 
driving frequency, cjl, the Hamiltonian reads: 



H 



hA c dla c 



+ :i Tr-Sz + hg{S+d c + S-al), 



(2) 



where A r 



wl and A q 



denote the de- 



tuning of the driving frequency, wl , from the cavity and 
qubit resonance frequencies, lo c and u> q , respectively, and 
g is the coupling strength between light and qubits. The 
cavity-field creation and annihilation operators are de- 
noted by, aj and a c , respectively, while Sk = Sj=i &k 

for k = +, — , z are sums of Pauli operators, ai r , for the 
two qubits. The coherent driving amplitude, /3, is nor- 
malized such that |/3 1 2 is the incident number of photons 
per second onto the input mirror, the field-decay rate 
of which is n\. Similarly, with re 2 being the field-decay 
rate of the exit mirror, the leakage of the cavity field 
is modeled by the decay operator, £1 = \/2rea c , in the 
decay part of Eq. (fT]), where re = rei + Population 
decay of qubit I and 2 with rate, 711 , can be modeled 

by £2 = y / 7j[<7l 1 \ £3 = % /7j[(rf 2 \ respectively, whereas 
collision-like phase decay of each dipole moment is mod- 
eled by £4 = —^=a^ and £5 = — }==ai 2 \ where r is the 
mean waiting time between the phase-disrupting events. 
The output field is given by 14] ot — V2re2a c — v, where 
the vacuum field, v, reflected from the exit mirror pre- 
serves the operator commutation relations but gives no 
further contribution at zero temperature (methods for 
treating finite-temperature environments are outlined in 
Ref. |l5|). The balanced homodyne detection setup mixes 
the output field and the local oscillator field, «LO; lead- 
ing to the differential [l3j and summed photo-currents 
(in units of electrons per second) : 



i+(t) = J7|a L o| 2 , 



(3) 



where the field-quadrature operator, Xg = ^{d c e~ ie + 
ale l6 ), depends on the relative phase, 9, of the local os- 
cillator. The operator Xg is connected to the formal- 
ism of Eq. (TTJ) when d is defined by d = \/2n2d c e~ ie , 
i.e. d + <v = 2\/2K2Xg. By defining the normalized dif- 
ferential photo-current, I(t) = i-(t)/ y/i+(t), we obtain: 



I(t) = 2^2^j(X e (t)) +£(t) = Jrj(d + <fl)+t(t). (4) 

With the notation, (A, B) = (AB) - (A) (B) , the correla- 
tion function, R(t; r) — (I(t + r), 7(f)) e, of the normal- 
ized differential photo-current is given by [13j ]: 

R{t;r) = 8K 2 r}(:X(6,t + T),X(6,t):) +S(t), (5) 
where ":" means normal-ordering of the field operators. 



A. Adiabatic elimination of the cavity-field 
variables 

The above dynamical equations are very general and 
can be simplified in our case of k 3> g by adiabatically 
eliminating the cavity-field variables. Our elimination 
procedure varies only slightly from previous works (see 
e.g. [l6l ]) and hence only the main steps are given: The 
cavity- field operator is written as a c = a c + a' c , where 
a c = corresponds to the mean cavity-field in ab- 

sence of qubits. Next, transform the master equation 
to the frame rotating at the qubit resonance frequency, 
w q , and eliminate adiabatically a' c . The resulting master 
equation is then transformed back to the frame rotating 
at ojl with the effective qubit Hamiltonian given by: 

g eff = ^S z + hg(a c S + +a*J-)- HA f" S + S ~ , (6) 
2 k 2 + A£ q 

where A cq = A c — A q . The qubit-decay operators, c m 
(with m — 2,..., 5), in the master equation (JlJ are 
maintained while the cavity-leakage operator is replaced 
by the correlated qubit operator, c\ — > ^/jpS-, where 

7 P = Jz+£i ■ The measurement operator, d, is replaced 

by: 

d eS = V2^(a c - %9 ^- )e-* e , (7) 
^ k + iA cq J 

which in turn from Eq. ((4]) leads to the photo-current: 

/(*) - sin(e - e K ) + (s y ) cos(e ~ e R )] 

2y / 2Ki] c ff [Re{a c } cos# + Im{a c } s'm9] + £,(t), 

' (8) 

where t&n0 K = S± = ±(S X ± iS y ), and r) eS = 

is the total detection efficiency accounting also for the 
non-detected fraction, ^ , of photons leaking through the 
left-hand mirror in Fig. [T] The photo-current correlation 
function turns into: 

R(t; t) = j p r) e «[(S + (t + r), S_(t)) + (S+(t),S-(t + r)> 

-(:5_(t + r),5_(i):)e- 2 ^)] +( 5(r), 

(9) 

where the normal-ordering is transferred from a c , aj to 

When inserted into the master equation (JXJ) , the two 
left-most terms in the Hamiltonian (|6]) together with the 
qubit-decay terms given by £2 , . . . , C5 corresponds exactly 
to the semi-classical description of light-matter interac- 
tions. In addition, the presence of the cavity introduces 
a Stark-shift term (right-most term in Eq. ([B])) and an 
additional, correlated spontaneous decay process by the 
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FIG. 2. (Color online) (a) The photo-current (upper part) 
and the singlet-state overlap (lower part) as a function of 
time for two individual instances of the simulations. One 
(blue crosses) eventually collapses into | — ) whereas the other 
(black circles) collapses into the triplet space. Panel (b) shows 
for each integration time, T, the distribution of the singlet- 
state overlap on a relative scale set by the shaded bar. 

ci-operator (the so-called Purcell effect [ijj )■ The elimi- 
nation of the cavity field is a good approximation when- 
ever k 3> 5,X:7||: r_1 ? where \ = 2ga c is the resonant 
Rabi frequency of the qubits. For the remaining part 
of this manuscript we assume the adiabatic elimination 
of the cavity field to be in effect when referring to the 
master equation ([1]). 

This parameter regime is relevant for transmission- 
wave-guide resonators [18|-|20j , in which the coupling of N 
electronic spins has reached the strong-coupling regime, 
\/N g ^ k,7||,t _1 , while typical values for the coupling 
parameter, g/2ir, to a single spin could be extended to, 
say, 300 Hertz. As we shall learn in Sec. [V] the qubit- 
decoherence rate must be significantly smaller than 7 P , 
which could be realized by coupling, e.g., single atomic 
ions [21] to such wave guides. 

B. Optimal strategy for entanglement generation 

The stochastic master equation ([T]) establishes the con- 
nection between, on one side, our knowledge of the quan- 
tum state described mathematically by p, and on the 
other side, the measurement record, I(t). This connec- 
tion can be established in two different ways. (I) From a 
theoretical perspective, Eq. JT]) presents a tool for simu- 
lating the realizations of measurements. The stochastic 
function, <!;(£), is then generated by the computer soft- 
ware and gives rise to the particular instance of p(t) 
and eventually the photo-current given by Eq. ([8]). (II) 
From an experimental perspective, the stochastic mas- 
ter equation can be employed for analyzing real exper- 
iments in which the photo-current, I(t), has been mea- 
sured. The function £(i) then represents the random- 
ness of the measurement process, and by continuously 
updating the coupled equations ^ and ([8]), the density 
matrix, p(t), will always correspond to the best obtain- 
able knowledge of the two-qubit quantum state. Even 
though the present work employs method (I) for simulat- 
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ing the entanglement-generating process, the full access 
to p gives the possibility to judge how well method (II) 
would work in experiment. 

The evolution of p and I(t) during the measurement 
process is exemplified in Fig. [2j For now we focus on the 
qualitative features and we defer a discussion of the spe- 
cific physical parameters to our detailed presentation of 
results in Sec. IIVI In panel (a) two instances of the sim- 
ulations have been shown; one which collapses into |— ) 
(blue crosses), and one which does not (black circles). 
Despite the fact that the two photo-current examples are 
both quite noisy, they do contain enough information in 
order to increase the knowledge of the singlet-state over- 
lap, which eventually becomes zero or unity as shown 
in the lower part of panel (a). In panel (b) the singlet- 
state-overlap distribution is shown versus time (based on 
10,000 simulations). This overlap is initialized at 50% but 
soon attains a much broader distribution. However, after 
few times 7" 1 of measurement, the overlap-distribution 
bifurcates into sharp peaks at zero and unity. Hence, for 
a sufficiently long measurement time the continuous evo- 
lution of p effectively facilitates the desired wave-function 
collapse. 

The optimal entanglement-generating protocol simply 
uses p to check to which degree the singlet state, |— ), 
has been realized. By requiring a minimum value, -F m i n , 
for the state overlap, the acceptance criterion for a given 
quantum state then becomes (— ) > -F m i n . 



C. Practical strategy for entanglement generation 



The optimal strategy discussed above requires knowl- 
edge of all the physical parameters, Kx, K2, A c , A q , g, 
7|| , r, in addition to sufficient data processing capabil- 
ity. This is indeed possible but might be impractical in 
reality, and hence some more robust but less accurate 
procedures for establishing whether the singlet state has 
been prepared are desired. To this end we shall consider 
the two integrated, dimensionless measurement signals: 



c</; = 7t I " I[t)dtl 



Clockin — m I 

Tti Jo 



(10) 



dt. (ii) 



In experiment, the former of these corresponds to a 
simple integration of the photo-current, while the lat- 
ter corresponds to inserting the photo-current signal 
into a lock-in amplifier with demodulation frequency, 
A;, and time constant, t;, and integrating for the mea- 
surement time, T, the squared modulus-output value, 
R 2 = X 2 + Y 2 (X and Y are the measured in-phase and 
in-quadrature amplitudes of the signal at frequency, A;). 
Now, our task is to devise conditions, e.g., \C{i)\ < C*/) r 

or Clockin < Cioddn' ^° acce pt the quantum state as be- 
ing sufficiently close to |— ). When simulating the en- 
tire entanglement-generating process, the procedure is 



repeated iVtotai times, and if A^ accept of these simulation 
runs lead to acceptance of the quantum state, we define 
the success probability as p suc = Accept /ATtotai- At the 
same time, the fidelity F — — - (~ I P \~~) measures 
the average occupation of the spin-singlet state for the 
generated quantum states, where the sum runs over the 
accepted density matrices, p, simulated by the master 
equation (fTJ). In principle, it should be possible to reach 
F = 1 with p suc — | (since the initial state has a 50% 
overlap with |— )). However, in practice a finite measure- 
ment time, qubit decoherence, and a non-optimal extrac- 
tion of information from I(t) reduces the fidelity obtained 
at a given p snc - An acceptance criterion is well-chosen if 
both F and p SU c attain high values. 



IV. ENTANGLEMENT GENERATION IN 
ABSENCE OF QUBIT DECAY 



This section is devoted to the generation of the spin- 
singlet state in absence of qubit-population and qubit- 
phase decay as modeled by the decay operators ca, . . . , eg, 
i.e. we take 711 = 0, r = 00. This simplifies the intro- 
duction of all the detailed concepts in the measurement 
scheme and defines the limits imposed solely by the mea- 
surement setup and by the chosen acceptance criteria. 
The influence of qubit decay is discussed in Sec. [V] 

In the numerical simulations the stochastic part of 
Eq. (TTJ) is integrated by the Milstein formula [22|, and 
p is evolved using time steps, dt, being 10 -2 times the 
characteristic decay time or oscillation time of the physi- 
cal variables. By repeating the simulations 10,000 times, 
the statistical spread on fidelity estimates is of the order 
of one percent. 

Without loss of generality, the phase of the driving 
field, f3, can be chosen such that a c and % are real. We 
shall also take A c = 0, i.e. f3 must then be real. By 
taking kx ~ and k% = K — n\ ~ k (i.e. the entire cav- 
ity decay takes place through the right-hand mirror in 
Fig-Hi tne effective quantum efficiency, i] c s, corresponds 
to that of the detector, r] e s ps rj. In experiment such a 
mirror asymmetry is not necessarily realistic, but adding 
a homodyne detection setup to the left-hand mirror out- 
put and combining the knowledge from all measurements 
would re-establish r7eff 77, and hence the choice Kx ~ 
just simplifies the simulations while maintaining the ex- 
perimental realism. The narrow qubit linewidth calls for 
A q = A cq <C k, and the correlated decay rate becomes 

2a 2 

7 P ss In the remaining part of this manuscript, all 
rates are measured relative to 7 p and time is measured 
in units of 7^ 1 . For the simulations we take specifically 

ensures the 



k — 50007 P and g = 50 7 P , such that 
validity of the adiabatic elimination. 
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FIG. 3. (Color online) All panels correspond to A q = 107 P , 
X = 1.65A q , 6 = — 7r/2, and 77 = 1. (a) Two instances of 
the integrated measurement signal, £m, as a function of the 
integration time, T. These examples are based on the photo- 
currents shown in Fig.[2ja), which led to the singlet state (blue 
crosses) or the triplet state (black circles). The red-dotted 
lines mark a chosen threshold condition, |C(i)| < C(/> r = 1-96, 
separating accepted and rejected states after the integration 
time, 7 P T = 10. Panel (b) shows from 10,000 simulations 
the distribution of the measurement signal, £/n , in case of 
HpH > 0.8 (crosses) and (-|p|-) < 0.2 (circles). The 
solid lines show Gaussian distributions with unit variance 
and with mean values given by the approximate estimate of 
Eq. (|12l) . The acceptance window from panel (a) is marked 
with red-dotted lines, (c) The obtained fidelity, F (dashed 
line), and success probability, p aU c (solid line), as a function 
of the threshold value, Ctj)i of the acceptance window. The 
red-dotted lines correspond to the choice of panels (a,b) lead- 
ing to p suc = 50% and F = 94%. (d) Solid line: The fidelity, 
F, as a function of the success probability, p suc . The red- 
dotted line marks p SU c = 50%. Blue-dashed line: Theoret- 
ical limit for completely separate singlet- and triplet-space- 
measurement signals and perfect state overlap. 



Measurement schemes 



DC-analysis of the photo-current 



Following the strategy presented in Sec. IIII C[ we con- 
sider first the use of the photo-current mean value signal 
from Eq. (jTUJ) for distinguishing between the singlet state 
and the triplet space. According to Eq. (|SJ) it is feasible 
to choose the local oscillator phase, 9 = — such that 

I(t) = \/jpVcs(S x ) does not contain a background 

contribution from the cavity field (note 9 K ps 0). The 



ensemble average of Eq. (ITU)) then becomes: 

(C(/))k 



7 P ??cff 



T 



((S x (t))) E dt 



^"f p Tr) e g{S x }ss (triplet) 
(singlet) 



(12) 



where the triplet-space steady-state value of S x is pre- 
dicted to be (details given in Appendix |A| : 



{S x )ss = 



-2xA q (7 



4A2- 



2X 2 ) 



( 7 2 4-4A2)( 7 2+A2 +X 2 ) 



3 Y 4' 



(13) 



and the approximation assumes that (S x (t)} corresponds 
to (S x )ss most of the time, i.e. j p T ^> 1. For the singlet 
state the 5-correlated nature of £(i) leads to the vari- 
ance, Var(£{/)) = 1, while the distribution is broader 

for the triplet space due to temporal variations in (S x ) 
(see Fig. ISJb) and the discussion below). Despite the 
crudeness of the approximation in Eq. (|12[) it is clear 
from the above discussion that an effective distinction 
between the singlet- and triplet-spaces is possible when 
JpTr) e g(S x )s S ^> 1, and that the parameters, A q and Xi 
should be optimized in order to maximize (S x )ss accord- 
ing to Eq. (fl"3)) . When A q 3> 7 P , the maximum value of 
Eq. (liU)) is (S^ss ps 1.52 obtained when x ~ 1.65A q (see 
also Fig. HJc)). 

Now, consider Fig. [3] exemplifying the entanglement- 
generating process. In panel (a) the integrated photo- 
current, £/n j of Eq. (ITU1) has been plotted for the two 
individual simulation runs, which were already discussed 
in Fig. [2ja). As time evolves, the black-circled curve 
shows an increasing value of | C<z> I ; which reflects the fact 
that the photo-current, I(t) (black circles), in the upper 
part of Fig. [Ha) has a mean value slightly below zero as 
a consequence of (S x ) being non-zero for a triplet-state 
simulation instance. In contrast, in Fig. [5Ja) the photo- 
current (blue crosses) representing a simulation instance 
ending up in the singlet state is closer to zero on aver- 
age, which again is reflected in the measurement signal 
(blue crosses) in Fig. [3Ja). Now, the practical accep- 
tance criterion consists simply of keeping a given state 
p provided that ends up at time 7 P T = 10 between 
the red-dotted lines in Fig. 02[a), i.e. if |C(/}| < C(/) r f° r a 
pre-selected value of Ca- 
using the entire set of simulations, the distribution of 
C(i) has been plotted in panel (b) showing a clear double- 
peak structure. By distinguishing between high (> 0.8) 
and low (< 0.2) overlap with |— ) , we clearly see that each 
peak corresponds to either the singlet or triplet space. 
The solid lines are Gaussian functions with unit variance 
(the shot-noise level of the homodyne-detection proce- 
dure) and mean values predicted by the crude approxi- 
mation of Eq. (fl~2"|) . The singlet-state (crosses) is modeled 
accurately since (S x ) = is exact and the only variation 
arises from the random shot noise of the measurement. 
However, for the triplet state the simulated distribution 
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is evidently broader and asymmetric — the additional 
width arises from the qubit dynamic evolution within the 
triplet-state manifold leading to a variation in S x . Again, 
the red-dotted lines depict the acceptance window, which 
clearly selects most of the singlet-state events; however, 
a small fraction of the undesired triplet-state occurrences 
are also included. This effect illustrates the fact that the 
experimentally simple acceptance criterion, |C(/)I < Ql)> 
is less accurate than the complete calculation of p dis- 
cussed in Sec. MI Bl In fact, if the optimal method of 
accepting states with (— \p\— ) > F min for some selected 
value of Fmin is used the undesired instances from the 
triplet space with poor singlet-state overlap would sim- 
ply not occur. 

The value of p suc and F can be calculated as a function 
of (*^ r , i.e. for various widths of the acceptance window, 
as shown in panel (c). Clearly, for small, increasing val- 
ues of C/ n r the success probability grows quickly without 
much degradation in fidelity since the acceptance win- 
dow selects predominantly the states with a high singlet- 
state overlap. When the 50-percent success probability is 
reached, a further increase of must incorporate some 
triplet-state instances with a loss of fidelity as a result. 
These observations can also be shown as an _F-versus- 
p suc plot, see panel (d). Here the ultimate limit (shown 
by a blue-dashed line) can be obtained if the singlet- and 
triplet-spaces present distribution functions like those of 
panel (b) but being entirely separate. 

The fidelity obtained at 50% success probability is 
shown in Fig. SJa) for various values of A q and \. The 
variation in this fidelity can then be compared to the 
triplet-state steady-state mean value, (S x )ss, which has 
been plotted in Fig. @|c) f° r the same parameter set- 
tings of A q and \- The correlation between these figures 
is evident, which confirms the simple picture discussed 
around Eq. (fT2"|) that the triplet-state imprint onto the 
photo-current must be maximized for optimizing the per- 
formance of the protocol. 

2. AC-analysis of the photo- current 

The spectrum of the photo-current is defined as: 

Sj(A) = — R(t;ry AT dT, (14) 

which generally depends on the time, t, in the initial 
transient regime but is time-independent in steady state. 
For the singlet state the spectrum is flat, Si (A) = w-, 
corresponding to the shot noise level of the homodyne 
detection apparatus. The triplet space is distinguished 
from the singlet state by identifying a spectral peak in 
this flat background. In similarity with the expectations 
discussed around Eq. (fl~2j) for the DC-analysis, we shall 
here use the steady-state spectrum of the triplet space to 
predict the optimization of A q , x, A;, and Tj for best 
performance of the measurement signal, Ciockin- As a 




FIG. 4. (Color online) Panels (a) and (b) show as a func- 
tion of the resonant Rabi frequency, \, the obtained fidelity, 
F, with a success probability, p suc = 50%. Panel (a): DC- 
analysis of I(t) with 9 = — 7r/2, and A q /7 P given by 10 (blue 
crosses), 1 (green circles), and 0.3 (red squares). Panel (b): 
AC-analysis of I(t) with 8 = 0, and A q /7 P equal to (blue 
crosses), 1 (green circles), 3 (red squared), and 10 (black di- 
amonds). The simulations corresponding to the blue crosses 
have been performed twice in order to depict the statisti- 
cal uncertainty of F. (c) The magnitude of the triplet-space 
steady-state value of S x as given by Eq. (|13[) for A q /7 P = 10 
(blue-solid line), 1 (green, dash-dotted line), and 0.3 (red- 
dashed line), (d) The height of the spectral peak in Si (A) 
when Aq/7 P = (blue-solid line), 1 (green-dash-dotted line), 
3 (red-dashed line), and 10 (black-dotted line). 

first step, consider the dynamical mean-value equations 
for S = [S x , S y , S Z ) T , which follow immediately from 
Eqs. (|AT|l - (|A"g|) with x being real: 

M = Qx <S)+ 7p [...], (15) 

where the coherent driving vector is given by Q = 
[x, 0,A q ] T , and the 7 p -term (with quadratic S- 
componcnts left out for clarity) tends to drive (S) toward 
the vector [0,0,-2], i.e., the state |gg). The modulation 
of the photo-current, through (S x ) and (S y ) according to 
Eq. (JSj) , is largest when the spin vector is allowed to sweep 
across the full sphere, i.e. we expect that x 3> 7p> A q is a 
good choice in order to maintain a significant level of ex- 
citation. We note from the expression of Q that (S) will 
primarily be spinning around the rc-axis, which leads to 
significant oscillations in (S y ). For this reason, the local- 
oscillator-phase choice, = 0, is natural. The oscillations 
in I(t) are then superposed on a constant background 
level oc Re{a c }. 

The spectrum of the simulated current, Ij — I(tj) 
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FIG. 5. (Color online) All graphs are based on 10,000 simu- 
lations with x = 107p> A q — 0, 7 P T = 10. (a) The simulated 
spectrum (dots) using the periodogram, Pi (A) [Eq. (|16[) ] on 
simulation instances with (— |/S|— ) < 0.2, compared to the an- 
alytical spectrum, Si (A) (solid line), of Eq. (|14[) . The dashed 
line at Ao = 9.897 P marks the maximum of S/(A), and V 
denotes the FWHM of Si (A), (b) The distribution of Ciockin 
distinguished by a high (> 0.8, crosses) and a low (< 0.2, cir- 
cles) singlet-state overlap. The lockin parameters, A; = Ao 
and ti — r , follow the characteristics of Si (A) from panel 
(a), and the red-dotted line marks an acceptance criterion, 
Ciockin < ClocMn) witn Psuc = 50%. Panels (c) and (d) show 
the obtained fidelity (with p suc = 50%) for varying lockin pa- 
rameters, A; and n, respectively. The dashed lines mark the 
parameters used in panel (b). Panels (e) and (f) show the 
variations in fidelity and success probability for varying ac- 
ceptance thresholds. The red-dotted lines mark the obtained 
fidelity F = 93% at p suc = 50%. 



where j runs over the discrete times separated by dt, 
can be conveniently estimated by the periodogram 



-Pr(A) 



dt 
2irn 



iAti 



(16) 



which is essentially the modulus square of the discrete 
Fourier transform of I(t). The front factor ensures the 



correct value, .Pr(A) = for the shot noise background, 
and we subtract from I(t) the constant contribution of 
the bare cavity ex Re{a c } prior to insertion into -Pj(A). 

Turning to the simulation, Fig. [HJa) shows the sim- 
ulated spectrum for the subset of instances, which col- 
lapse into the triplet space, in comparison to the expec- 
tation from Eq. (fl4|) . which in steady state is given by 
Eq. (|A8|) . Since the integration time, 7 p T = 10, is sig- 
nificantly larger than unity, this steady-state expression 
does in fact match the simulated curve very well. In the 
limit, x 3> A q ,7 p , numerical inspection of Eq. (| A8|) re- 
veals a Lorentzian peak centered around the generalized 

Rabi frequency, = 



X 



and with full-width at 



half maximum (FWHM) T = §7 P (for the solid curve 
in Fig. [S^a), the maximum is placed at A = 0.989J7 
with r = 1.547 P ). In order to distinguish between the 
singlet- and triplet-space part of the initial state, we 
must establish the absence or presence of this spectral 
peak in each individual simulation run. To this end, the 
lockin parameters of Eq. (ITT1) are chosen as A/ = Ao and 
ti = r _1 , and as can be seen from Fig. [5jb), the mea- 
surement signal, Ciockin, is indeed capable of separating 
the singlet state from the triplet space. The robustness 
of this procedure to errors in the lockin parameters is 
depicted in Fig. [Hc,d), which show that A; must obvi- 
ously match the position of the spectral peak with an 
accuracy set by T and that r/ w provides the best 
match to the bandwidth of the signal peak. The obtained 
fidelity and success probability while varying the accep- 
tance Criterion, Ciockin < Ciockin' can be seen in Fi g-[2I e :f)- 
These graphs are quite similar to the corresponding re- 
sults for the DC-measurements in Fig.[3Jc,d); however, if 
a low success probability is accepted, the obtained fidelity 
seems to be better. 

The optimization of the qubit-driving parameters, x 
and A q , are examined in Fig.[4jb,d). Panel (b) shows the 
obtained fidelity in various simulations runs, and there 
is a clear correlation with the calculated spectral-peak 
height of <S7(A) shown in panel (d). In similarity with 
the DC-analysis in Sec. IIV A 11 the present AC-analysis 
of I(t) is optimized in terms of x an d A q simply by 
maximizing the triplet-state steady-state spectral peak 
height, and Fig. |U[d) presents the practical condition, 
X ^ 10 • max(7 p , A q ), for this optimization. 



B. Practical versus optimal extraction of 
information 



In Fig. [6] the performance of the DC-analysis of I(t) 
as described in Sec. IIV A II (solid triangles) can be di- 
rectly compared to an optimal extraction of information 
from the full I(t) (solid circles) given the local-oscillator 
phase, 9 = —ir/2. Likewise, the lockin-based AC-analysis 
of I(t) [Sec. IIV A 2| shown with open triangles can be re- 
lated to an optimal information extraction (open circles) 
given the phase choice, 9 = 0. In comparison to the 
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FIG. 6. The obtained fidelity versus integration time when 
Psuc = 50% (panel a) or p suc = 10% (panel b). Fidelities 
obtained by optimal analysis based on (— ) > F m i n are 
shown for 6 = (open circles) and 6 = —iv/2 (solid circles). 
The lockin analysis, Ciockin < Ciockin f rom Sec. IIV A 21 gives 
rise to the open triangles, while the integrated-photo-current 
analysis, \C(i)\ < Cf[) from Sec. IIV A 11 leads to the solid 
triangles. The solid squares represent the displaced, time- 
weighted measurement signal described in the text. 



optimal extraction of information, the simple and more 
robust approaches require approximately twice the time 
for obtaining a given fidelity with a given success rate. 

For the DC-analysis protocol the detection record dur- 
ing the initial transient dynamics of duration ss 7" 1 does 
not bear much information since neither the singlet- nor 
triplet-space part of gives rise to a non-zero value 
of (S) in the initial time range, < t < Jp 1 ? m which 
the two qubit spins are pointing in opposite directions. 
Hence, by weighting the integral in Eq. (fT0|) by the func- 
tion (1 — e _7pt ), we do not loose information but a smaller 
amount of shot noise is accumulated in this transient part 
of the protocol. Combining this weighting procedure with 
an acceptance window shifted by 0.5 toward the right in 
Fig-GHb), we obtain the improved fidelities shown by solid 
squares in Fig. [5J and the DC-analysis protocol narrows 
in on the full calculation of p. 

The AC-analysis protocol is based on the correlation 
function Q, which contains quadratic moments of S and 
hence is able to deliver an oscillatory signal starting al- 
ready from t — 0. Considering Fig. [6l we ascribe this fact 
as the reason for the slightly better performance of spin- 
precession-based protocols (open symbols) in comparison 
to the spin-mean- value-based protocols (solid symbols) at 
short integration times. 



V. THE INFLUENCE OF DECOHERENCE 
PROCESSES 

This section estimates the effect of decoherence pro- 
cesses on the obtainable fidelity. We note that if such 
processes are strong, the optimum parameter settings as 
exemplified by Fig. 0] might change. Instead of perform- 
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FIG. 7. (Color online) All panels show simulation results for 
A q = 107 P , x — 16.57 P , and 9 = —n/2. (a) The obtained 
fidelity when p suc = 50% using |C(i)| < C(i> while varying 
7||/7p = 1 x 10~ 3 (Red squares), 2 x 10 -3 (green circles), 
5 x 10~ 3 (blue crosses), 1 x 10~ 2 (black diamonds), 2 x 10~ 2 
(cyan tip-up triangles), 5 x 10~ 2 (magenta tip-down trian- 
gles), (b) The optimum fidelity when p suc = 50% for various 
characteristic decoherence times t c . Red squares: Popula- 
tion decay with r c = 7jT , green circles: Phase decoherence 
with t c = r, blue triangles: Inhomogeneous coupling strength 
with r c = Sx, black diamonds: Inhomogeneous qubit fre- 
quency with t c = Su! q . Open and closed symbols are obtained 



using < C/> r and (— ) > Fmin, respectively (this 



holds also for panels c and d). (c) Fidelity versus 7,, for 
p SU c = 50% (squares), p suc = 30% (circles), and p suc — 10% 
(diamonds) . Solid and dashed lines denote F opt = 1— ^Ayy/Tp 
and F op t = 1 — 711 /7p, respectively, (d) Fidelity versus #X _1 
for p suc = 50% (triangles) and p suc = 10% (squares). Solid 
and dashed lines correspond to slopes of —1 and —2, re- 
spectively. In panels (e) and (f) \ip) = | — ) at t = and 
7||/7p = 0.01. (e) singlet-state overlap at 7 P T = 1 (blue 
crosses), 2 (green circles), 5 (red squares), and 20 (black dia- 
monds), (f) average value of (— ) for all incidences (blue 
crosses) and for those with (— \p\— } > 0.5 (green circles). Red 
squares: Most probable value of (— \p\ — ). Black diamonds: 
Fraction of states with (— \p\ — ) > 0.5. Solid line: exp(— 7||T). 
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ing a full-scale analysis of such possible changes, we sim- 
ply add decoherence processes but keep the measurement 
protocols and acceptance criteria. The analysis presents 
a lower bound of the obtainable fidelities an must be a 
good approximation in the limit of high fidelities. Only 
the local-oscillator-phase choice of 9 = —ir/2 relevant for 
the DC-analysis protocol is discussed below — the case 
of 8 = presents similar features. 

To exemplify the simulation procedure, a qubit- 
population decay is introduced in Fig. Eta) with varying 
values of 711 (modeled by the decay operators £2 and £3 
in Sec. IIII[) . Since the decay process deteriorates the de- 
sired singlet state in the long-integration-time limit there 
exists an optimum integration time and a correspond- 
ing optimum fidelity. This fidelity has been plotted in 
Fig. [3(b) — see the figure caption for simulation details. 
In a similar manner, the effect of qubit phase decay can 
be modeled by assigning a finite value to r in the opera- 
tors £4 and £5. Furthermore, one may consider the case 
that the two qubits are not coupled in the exact same 
way to the cavity. A small inhomogeneity in the cou- 
pling strength for each qubit, g\ = g + 4^ and g^ — g— 4f, 
will lead to a slight difference in resonant Rabi frequency, 
Sx = 2a c (5(7, or alternatively, a small difference in qubit 
resonance frequency, <5w q , could be present. The effect 
of these non-ideal scenarios are compared in Fig. [71(b) 
showing that qubit population decay (red squares) and 
dephasing (green circles) behave in approximately the 
same way, while the inhomogeneities in Rabi frequency 
(blue triangles) or qubit detuning (black diamonds) fol- 
low their own distinct trend. These observations can be 
complemented by the equations of motion for the singlet- 
state population (only the deterministic part, i.e. take 

v = oy. 



dp— t - 
dt 



1 A 1 
- +711 P-- + -P+-+ +7||Po 



T 

iSuj. 



2 

iSx 



Pss- 



p-,+) 

- Pee - 



(17) 



P-, 



V2 



V2 



Evidently, the two rates 711 and r _1 enter on the same 
footing and are responsible for the de-population of the 
singlet state. In the high-fidelity limit {p~- ~ 1, 
P+.+ 7 Pcc,cc « 0) one would expect the infidelity, 1 - F, 
to increase linearly with these rates. Likewise, the in- 
homogeneities parametrized by Suj q and Sx seem to be 
comparable in effect — they attempt to drive coherently 
the population from the singlet state toward |+) (the 
(5o; q -term) or ^(|gg) — |ee)) (the Jx-term). In the high- 
fidelity limit the coherence terms, Pgg,- and p cc .-, 
must be polarized before the singlet-state population can 
be driven, and hence the infidelity is expected to increase 
quadratically with Sx or <5w q . In order to exemplify these 
scaling behaviors, consider the red squares and blue tri- 
angles of Fig. EJb), which have been re-plotted in pan- 
els (c) and (d), respectively. These data scale roughly 



as the solid lines, which in the double-logarithmic plots 
have slopes — \ and —1 in panels (c) and (d), respec- 
tively. This does not correspond to the scaling behavior 
discussed above; however, by accepting a smaller success 
probability, the fidelity increases and the expected scal- 
ing is found in the high-fidelity limit. In fact, for the 
best case shown in panel (c) with optimal information 
extraction and p suc = 10% (solid red diamonds), the in- 
fidelity becomes, 1 — F opt ~ T||/Tp- At nrs t glance this is 
surprising since the integration time for obtaining high 
fidelities is typically exceeding 107" 1 as exemplified in 
panel (a); however, the infidelity-penalty is not exceed- 
ing 107||/7 P but is equal to roughly one unit of 711 /7 P - 
A thorough examination of the continuous measurement 
process in the high-fidelity limit is required to explain 
this observation: Consider panels (e,f) based on 10,000 
simulations in presence of qubit-population decay and us- 
ing the singlet-state as the initial state, i.e. (— \p\— ) = 1 
at t = 0. Panel (e) shows the distribution of (— \p\ — } 
for various times, which is seen to increase in width dur- 
ing the first few 7 p T but then settles to an almost con- 
stant distribution (compare red squares and black dia- 
monds). However, a sharp feature is emerging around 
(— \p\— ) = showing an increasing population within the 
triplet space. Now, the ensemble mean value of (— \p\— ) 
based on these distributions is calculated and plotted in 
panel (f) (blue crosses). The dynamics of this ensemble 
mean value is governed simply by the deterministic part 
of the master equation (|T|) (since (£(£))e = 0) which for 
the p__-component is given by Eq. (|17[) . In the high- 
fidelity limit the presence of population decay leads to 
3p g t '~ sa — 7||/7_,_, the solution of which is the solid line 
in panel (f) confirming the simulations. Nonetheless, the 
distribution is strongly peaked around zero and unity. 
In comparison, a drastically different sub-cnscmble mean 
value of (— \p\— ) can be obtained if we are able to select 
the best singlet-state candidates. For instance, by con- 
ditioning the sub-ensemble on (— \p\— ) > 0.5, the green 
circles are obtained in panel (f). Furthermore, the most 
probable value of (— \p\— ) (i.e. the maximum location 
of the distributions in panel (e)) turns out as the red 
squares in panel (f). Clearly, after a few times 7 p T these 
conditioned observables settle to a steady value. This 
does not contradict the fact that the state is decaying — 
the fraction of states with (— \p\— ) > 0.5 is decreasing 
steadily as shown by black diamonds in panel (f). To re- 
capitulate the above discussion: Once the singlet state 
has been established the continuous measurement either 
preserves it with a high fidelity or the state jumps into 
the triplet space. Since the characteristic time for updat- 
ing our knowledge is Jp 1 , the infidelity of the preserved 
singlet state is approximately 7||/7 P , which explains the 
steady value of the red squares in panel (f) and the ob- 
served dashed-line scaling in panel (c) (since exactly the 
very best states are selected in this case). A similar set of 
arguments can be made for the high-fidelity limit in panel 
(d) for the case of inhomogeneous cavity-qubit coupling, 
and we note in both cases that the dashed-line slope is 
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in general maintained also for the DC-analysis-based ac- 
ceptance criterion, \ < Cm- We note that by select- 
ing the best states (in particular when using the optimal 
extraction of information) the observed infidelity is de- 
termined by the ability of the measurement to preserve 
the state and not sensitive to the statistics of the finite 
number of simulations. 

In the cases of qubit-population or phase decay, the 
product 7 p t c is essentially equal to (up to factors of 

two) the cooperativity parameter C = where j± = 

— + In the limit of maximum success probabil- 

ity, p suc — 0.5, the above observations conclude that 
1 — F op t f=a which is similar to the figure of merit 

for deterministic protocols of cavity QED [23| . The mea- 
surement process seems to counteract effectively the loss 
of coherence inherent in the bad-cavity limit, and by ac- 
cepting a moderately lower success probability, the states 
with an infidelity of 1 — F opt ~ ^ can be conditionally 
prepared. An infidelity scaling as ^ was also obtained in 
the heralded protocol of Ref. [23| . 



VI. DISCUSSION 

From a fundamental perspective, the correlated decay 
with rate 7 P is the key mechanism for the entanglement- 
generation protocols. The original product state with op- 
posite qubit-spins gives rise to a 50% overlap with |— ), 
which allows the measurement process to induce a col- 
lapse into the desired singlet state with a high success 
probability. However, the information of the two-qubit 
state must be transferred to the cavity field and subse- 
quently leave the cavity before reaching the homodyne- 
detection apparatus, and it is exactly the decay rate, 7 P , 
which describes the combined rate of this information 
flow. When the cavity field is adiabatically eliminated 
the decay rate, 7 P , materializes explicitly as a strength 
parameter in the homodyne-detection photo-current, see 
e.g. Eqs. (0 and (JH). Hence, in our analysis it does not 
make sense to consider the particular case of 7 p = (as 
was done in Figs. 2 and 3 of Ref. [Hj]) since no informa- 
tion is gained by the measurement. 

Even though the measurement process is continuous, 
the discussion in Sec. [V] revealed an effective jump-like 
behavior of the quantum state, and we also remind that 
discreteness is regained in the long-integration-time limit 
as exemplified clearly by Fig. [2] These observations are 
not only of fundamental importance — the detection sig- 
nal effectively monitors any unwanted transitions into the 
triplet space. A simple feedback can thus be implemented 
in order to establish and maintain a high spin-singlet 
overlap in a continuous operating mode of the experi- 
ment. We imagine that a continuous feedback strategy 
along the lines of Ref. [l6[ can be developed taking into 
account the known imprint of the triplet-state compo- 
nents onto the photo-current. 

The entanglement-generation protocol relies heavily on 



the fact that the |— )-state does not couple to the cavity 
field and the interaction with the field effectively imple- 
ments a QND measurement [24| of the projection op- 
erator q = | — ) (— |. This apparent QND character of 
the protocol does not rely on the adiabatic elimination 
and the protocol should be applicable outside the bad- 
cavity limit. However, the optimization considerations 
of Sec. IIVI and the fidelity analysis of Sec. fVl require the 
adiabatic approximation to be valid. 



VII. CONCLUSION 

A measurement-based entanglement-generating proto- 
col has been established for two qubits residing in a bad 
cavity. The separation of an initial state into either the 
singlet or triplet space facilitates the establishment of 
entanglement, which is done optimally by a stochastic- 
master-equation approach. In addition, two practical 
methods have been discussed: (1) The use of the inte- 
grated photo-current mean value and (2) a lock-in-based 
identification of oscillations in the photo-current. The 
optimization of these methods for best performance can 
be understood simply as maximizing the imprint of the 
triplet-space part onto the photo-current. The influence 
of qubit dissipation and inhomogeneities has been an- 
alyzed such that obtainable fidelities can be estimated 
from relevant experimental parameters. 
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Appendix A: Steady-state properties of the triplet 
space 

In this appendix the photo-current, I(t), and its spec- 
trum, Si (A), is calculated under the assumption that 
the two qubits have reached steady state and also as- 
suming that qubit-decay processes are absent (7y = 0, 
t = 00). Under the latter assumption the singlet and 
triplet spaces are decoupled and we consider here only 
the three-dimensional triplet space. We shall also ignore 
the information gain from the photo-current (i.e. leav- 
ing out the measurement-super-operator part of Eq. ([Pjl) 
in order to establish an a-priori prediction of the photo- 
current. 

Considering the master equation with the choices made 
above, the dynamical equations for mean values of 5+, 
S-, and S z read (with \ rea l an d neglecting the last, 
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Stark-shift term of Eq. ©): 



(Al) 



^ = -iA q {S-) + l f(S z ) + ^{§J-) : (A2) 



= -i X [(S+) - (S-)} - 2 lp (S+S-), 



(A3) 



d(S z 
dt 



In Eq. (|A3|) the second line is valid in the special 
case of the triplet space since the two operators 0\ = 
S+S- and 6 2 = 2 + \S Z 



-> -i X [(S+) - (S-)]+7 P [-(S 2 z ) - (&) -4]- 



\S 2 act identically on the 



triplet-state basis set: Oj |ee) = 2|ee), Oj |+) = 2|+), 
Oj |gg) = 0, for 7 = 1,2. We note from all three 
equations above that the right-hand sides contain mean 



values of the quadratic operators, S + S z , S Z S-, and 
S 2 ., and to proceed the time-derivative of these mean 
values must be calculated. In turn, cubic operators 
are introduced and the set of equations seems endless. 
However, as exemplified above for Eq. (|A3j) . by em- 
ploying operator identities valid in particular for the 
triplet space, the equations become closed within an 
eight-dimensional space (corresponding to the number 
of free parameters in the triplet-space density matrix). 
Considering the column vector of mean-values: x = 
[{§+), (S-), (S z ), (S + S z ), (Sj-), (S 2 + ), (S 2 ), (S 2 )] T , the 
dynamical equations become after some algebra: 



<9x 



(A4) 



where the matrix, A, and the column vector, b, are given 
by: 
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The steady-state value of x is now simply given by xss = 
A _1 b, and since S x = S+ + 5_ the result of Eq. (|T~5|) 
follows after some algebra. 

Turning to the spectrum, Si (A), of the photo-current 
as defined in Eq. (fH|) , we note that correlation functions 
between S+ and S- are required according to Eq. ([9]). 
To this end, define first a new vector: 



(S + (t + r)S_(t)) 
(£_(< + r)£_(t)) 

(S z (t + T)S_(t)) 
(S+(t + T)S z (t + T)S-(t)) 
(S a (t + T)S-(t + T)§-(t)) 

r)S-(t)) 
r)S-(t)) 
r)S-(t)) . 



(Sl(t 

(s 2 -(t 
(s 2 z (t- 



(A6) 



by S-(t) inside the mean value brackets (•••). According 
to the quantum regression theorem, the time-evolution of 
y(t, t) for r > follows the exact same equation as x(t), 
- Ay(t, r)-b{S-(t)), which has the solution: 



dr 



i.e. this is the vector x evaluated at t + r and multiplied 



ty(i,r)-y(t,oo)] = e AT [y(t,0) - y(t,w)l (A7) 



where y(i,oo) = A- 1 b(5_(t)> = x ss (S*_(t)). The first 
two^entries of [y(t, r) — y(£,oo)] are equal to (S+(t + 
r),S-(t)) and (S-(t + r), S-(t)), respectively, which 
(together with their complex conjugates) are exactly 
the terms required in Eq. @. Since in steady state, 
R(t; — t) = R(t; r), the spectrum can be calculated con- 
veniently as an integral over r > 0: 
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(A8) 



with C, d, and v denned as 



{v T [(A + zA) 



(A-zA)- 1 ][y(t,0)-y(< ;W )]+c.c.} 

(A - iA)- x ][(C - (5_»A- 1 b + d] + c.c.} , 
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(A9) 



In the second line of Eq. (|A8[) the row-vector, v , ex- 
tracts the two first entries of [y(t, r) — y(t, oo)} and 
multiplies these by the appropriate weights according 
to Eq. ©. In the third line the solution (|A7[) is used 
and the r-integration of e Ar cos(Ar) is carried out. The 
last line includes (in similarity with the discussion below 
Eq. (|A3I) ) the re-expression of the quadratic and cubic 
terms of y by the linear and quadratic terms of x (in 
steady state) as a particular property of the triplet space: 
y(i,0) = Cx(t)+d. 

For completeness, we present a single-qubit version of 
the above results in the more general case that decoher- 
ence processes and the Stark-shift term of Eq. ([6]) are 
included. The vectors x, y, and v are restricted to the 
first three terms only, and the correlated-decay operator 
reduces to C\ — ./7p<7_ presenting effectively an extra 
qubit-population decay channel. The last line of Eq. (|A8j) 



I 

remains valid when replacing: 
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(A10) 



where 7^ 



A, 



- \ with 7 |f ff = 7|| + 7p , and Af = 
The single-qubit version of Eq. (|13|) becomes: 



xAf 



(Af) 2 + ( 7 f) 2 (l + 



71" 7 
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